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We study the meta-stable states in high temperature phase of QCD characterised by nonzero 

expectation values for the imaginary part of the Polyakov loop. We consider Nf = 2, 3 dynamical 

00 staggered quarks, and carry out simulations at various values of the coupling f3 to observe these 

i— H 

states. In particular, we find the value of the coupling (/3 m ) above which the meta-stable states 

-i— > 

appear. The resulting value of /3 m corresponds to temperature T m > 750 MeV for Nf = 2. 
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I. INTRODUCTION 



Experimental and theoretical studies of matter under extreme conditions is one of the 
active fields of research in recent times. Such a state of matter is created in ultra-relativistic 
heavy-ion collisions. The fireball created in these collisions leads to a state of deconfined 
quarks and gluons (Quark-gluon plasma). With the increase in the collision energy, the fire- 
ball not only crosses the confinement-deconfinement transition temperature T c but probes 
deeper into the deconfined phase. For example, in the heavy-ion collisions at LHC the initial 
fireball temperature is expected to go up to 5T C [1]. The increase in fireball temperature 
will result in the observation of new signals bearing the properties of the system at higher 
temperatures. Hence, it is important to study any possible prominent changes in the prop- 
erties of the medium in the deconfinement phase away from T c , which can be observed in 
experiments. In this context, we plan to study the explicit breaking of Z{3) symmetry, and 
the meta-stable states associated with it. 

In the pure £77(3) gauge theory, the deconfined phase exists in three degenerate states 
characterised by three different values of the Polyakov loop. These three states are related 
to each other via the Z{3) rotations. So, in the deconfined phase, Z{3) symmetry is spon- 
taneously broken. In the confinement phase, the Polyakov loop average vanishes restoring 
the Z(3) symmetry. For QCD with dynamical fermions, the Z(3) symmetry is explicitly 
broken, and the degeneracy between the three states is lifted [2]. Only the state for which 
the expectation value of Polyakov loop is real becomes the ground state. It is not clear what 
happens to the other two states. For asymptotically large temperatures, one expects the 
gluons to dominate so that the effects of quarks can at most make the other two states (with 
Polyakov loop phase angle ±27r/3) meta-stable. 

We would like to emphasize that the Z(3) meta-stable states are not the results of meta- 
stability near any first order transition. They are different from the meta-stable states that 
one observes near a first order phase transition. A closer analogy of the Z{3) meta-stable 
states would be a state of magnetization anti-parallel to the external field. 

The Z{3) meta-stable states are expected to play important role both in the context of 
heavy-ion collision, and in the early Universe. If these states indeed exist just above T c , then 
they can have significant effects on the medium properties [1]. If a fluctuation in the form of 
a meta-stable bubble in the back-ground of the stable phase costs free energy of the order of 
the temperature scale T, then such fluctuations will be present in the system. In the heavy- 
ion collision it is possible that the whole fireball may thermalise to one of these meta-stable 
states. This state will then decay through a first order phase transition even before the 
system cools down below T c . It is also possible that the meta-stable phases of super-horizon 
size may occur in the early Universe, and decay through bubble nucleation [3, 4]. 

There are several studies on Z(3) meta-stable states at high temperatures in the presence 
of quarks. It has been shown that the contribution of massless quarks to the one loop 
effective potential leads to meta-stable states for T > T c [2, 5]. These states have also 
been observed above the deconfinement transition in the Nambu-Jona-Lasinio model [6]. 
There are only a very few lattice QCD studies on these meta-stable states. A lattice QCD 
study with fermions in the sextet representation has found meta-stable states, characterised 
by phase angle (±27r/3, 7r), close to T c in the deconfinement phase [7]. However, in this 
work, we have considered Nf = 2,3 staggered fermions in the fundamental representation 
to look for the meta-stable states. Contrary to the previous studies, we find that the meta- 
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stable states do not exist in the neighbourhood of T c , but for temperatures T m > 750 MeV. 
Though this temperature may not be reached at the SPS and RHIC experiments, there 
is a possibility that these states are accessible at LHC. The fact that T m > T c may have 
important consequence (s) for the early Universe. 

The paper is organized as follows. In section II, we review the Z{3) symmetry and meta- 
stable states in the presence of dynamical fermions. Our lattice simulation techniques are 
discussed in section III. Results are presented in section IV. We present our conclusions in 
section V. 



II. Z(3) SYMMETRY IN THE PRESENCE OF DYNAMICAL QUARKS AND 
META-STABLE STATES 



In this section, we briefly discuss the pure SU(N) gauge theory at finite temperature, 
and their symmetries. Later, we focus on the case with dynamical quarks, and discuss their 
effects on the Z(N) symmetry [8, 9]. 

We start with the pure SU{N) gauge theory at finite temperature, T = In this 
case, one uses static fundamental charges (infinitely heavy test quarks) to probe into the 
dynamics of the pure glue system. The static fundamental charges are described by the 
Polyakov loop [10, 11] which is defined as the trace of the thermal Wilson line, 



L(f) = lTrW(x), 



(1) 



with the thermal Wilson line operator, L(x), defined as, 



W(x) = P exp 



ig / A (x,r)dr 
Jo 



(2) 



The expectation value of Polyakov loop (L(x)) is an order parameter for confinement- 
deconfinement transition in the pure glue theory. Here, P denotes path ordering of the 
exponential, g is the gauge coupling, and /3 T — 1/T denotes the extent of Euclidean time. 
A^ l {x ) t) = A®(x, r)X a is the vector potential. The A a are the TV 2 — 1 Hermitian generators of 
the SU (N) algebra in the fundamental representation. Aq(x : t) is the time component of the 
vector potential at spatial position x and Euclidean time r. The gauge fields A^{x^r) obey 
periodic boundary conditions in the Euclidean time direction, Ap{x, /3 T ) = A^x^O). These 
boundary conditions are maintained by a group of non-trivial gauge transformations [8] that 
are periodic up to a constant twist matrix, z G SU(N), 

g(x,f3 T ) = zg(x,0). (3) 



These matrices z form the center Z(N) of the gauge group SU (iV), where Z(N) is a cyclic 
group of order N . Thus the pure SU (N) gauge theory at finite temperature has the complete 
symmetry Q x Z{N), where Q is the group of strictly periodic gauge transformations. 

However, the Polyakov loop (which characterises the phases of pure gauge theory) trans- 
forms non-trivially under the Z(N) transformations, though the Euclidean action is invariant 
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under the transformations. Under the global Z(N) symmetry transformations, L(x) trans- 
forms as 

L(x)^zL(x), (4) 



where 



z = exp(2mn/N)l G Z(N), n G {0,1,2,- •• 



TV - 1}. 



(5) 



For temperatures above the critical temperature T c , the high temperature phase or the 
deconfining phase is characterised by (L(x)) ^ 0, corresponding to the finite free energy 
of an isolated heavy test quark, and thus breaks the Z{N) symmetry spontaneously. At 
temperatures below T c (in the confining phase), (L(x)) = 0, thereby restoring the Z{N) 
symmetry [9]. 

We now come to the effect of matter fields in fundamental representation on Z(N) sym- 
metry. In the fundamental representation of SU{N) ) the quark fields transform as, 



As a result the fermion part of the action is not invariant under the twisted transformations 
pertaining to the Z(N) symmetry. Thus in full QCD, the fermion part of the action breaks 
the Z(N) symmetry while the pure gauge part respects it. 

The effect of quarks on the Z{N) symmetry has been discussed in detail in [5, 12]. It 
has been advocated that one can take the effect of quarks in terms of explicit breaking of 
Z(3) symmetry [12-14]. This leads to a unique ground state and two (possible) meta-stable 
states in the deconfinement phase [2]. It has also been argued that the effects of quarks in 
terms of explicit symmetry breaking may be small, and the pure glue theory may be a good 
approximation [14]. However, we will see later that the explicit symmetry breaking due to 
quarks is too large for the meta-stable states to exist in the neighbourhood of T c . 



III. NUMERICAL PARAMETERS AND STUDIES 

In principle meta-stable states should ocuur in any simulation run. Given enough time 
the system will explore all of the configuration space. But the probability that an arbitrary 
initial configuration temporarily "thermalising" to a meta-stable state is very small. Once 
the system thermalises to the absolute ground state the probability that it fluctuates to 
the meta-stable state decreases exponentialy with the system size. So it is practical to 
choose initial configurations carefully so that they thermalise to the meta-stable states before 
decaying to the ground state. In order to have the system thermalized to a meta-stable state 
(ReL < and |ImL| > 0), the trial configuration should not be far away from this state. 
Otherwise, it will thermalize to the absolute ground state. Since a random configuration 
(termed as "fresh" in MILC code) results in a stable state (ReL > 0), we need to select the 
initial configuration appropriately. As the meta-stable states arise from the Z{3) symmetry 
so they are expected to be close to the Z(3) rotation of the absolute ground state. This 
fact can be used to find an initial configuration which may thermalize to the meta-stable 
state. To generate the gauge field configurations, we have used the MILC code which uses 
the standard Hybrid R algorithm [15]. We have used the similar simulation parameters as 



(6) 



in [7, 16]. 
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As a first step, we have performed a pure gauge calculation on a lattice of size 16 3 x 4 
near critical temperature (/? = f3 c = 5.6925 [17]) starting with a "fresh lattice". We then 
have selected one configuration each for the the following cases- 

(i) ReL < and ImL > 0, i.e. 9 ~ 2tt/3. 

(ii) ReL < and ImL < 0, i.e. 9 ~ -2tt/3. 
(hi) ReL <C and ImL ^ 0, i.e. 9 ~ tt. 

(iv) ReL > and ImL - 0, i.e. = 0. 

The case (iii) has been taken into consideration additionally to study whether a meta-stable 
state could exist at 9 = tv [7] , and the case (iv) to check whether we obtain the similar stable 
state if we start with a "fresh lattice" . Another method, we have employed, is by using an 

initial configuration with all the temporal links on a fixed time slice set to e-^^/^I. The 
pure gauge calculations have been performed using the MILC code. We have considered 4 
over-relaxation steps for each of the heat-bath iteration for updating a gauge configuration. 

We have used each of the gauge configurations (one for each 9) so obtained as an initial 
configuration to thermalize, and calculate the Polyakov loop in presence of 2 and 3-flavor 
quarks for a series of (3 values, 5.2 < f3 < 6.0. For this purpose, we have again used the 
MILC code with dynamical staggered fermion action. The numerical calculations have been 
performed with quark mass am u ^ = 0.01, so that m u ^/T = 0.04. Our micro-canonical time 
step size has been At ~ 0.01, and the trajectory length (micro-canonical time step x steps 
per trajectory) has been r ~ 0.8. Each gauge configuration has been analysed after 10 
heat-bath iterations. For each /3, we have collected 2500 gauge configurations. 



IV. RESULTS AND DISCUSSIONS 

A. How to observe the meta-stable states? 

In this study the Polyakov loop expectation value (L) has been used to charecterise 
different possible states of the system. In the absolute ground state (L) is real and positive. 

2n 

In the meta-stable states (L) is complex, with phase of (L) close to ±— . 

In a typical Monte-Carlo simulation, a sequence of statistically significant configurations 
is generated. This sequence constitutes the Monte-Carlo history. In order to reduce the auto- 
correlations between the consecutive configurations, measurements are perfomed on every 
10th configuration in the Monte-Carlo history. In all the different simulations, the initial 
configuration is thermalised in a few hundred Monte-Carlo steps. During the thermalisation, 
all the measured observables change almost monotonically. Once the system is thermalised, 
all the physical observables fluctuate around their respective average values. We compute 
the histogram (probability distribution) of the Polyakov loop, P(L). The meta-stable states 
are local minima in the free energy which are stable against small fluctuations. As a result, 
they appear as well defined peaks in P(L). Distinct Multiple peaks in P(L) correspond to 
different possible states of the system. 
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FIG. 1. Histograms of ImL vs. no. of trajectories (a) Single peak at ImL 7^ indicating a meta 
stable state only, (b) Double peak at ImL ^ and ImL = indicationg a meta stable state 
decaying into the ground state. 



In Fig. 1(a), we show the histogram of ImL (imaginary part of L) for two simulations 
with different random number sequences. The two P(ImL) distributions are similar given 
the small length (a few thousands) of the Monte-Carlo histories. The peak positions are at a 
non-zero value of ImL ~ —0.65 which imply a meta-stable state. This being a meta-stable 
state, it will decay to the ground state if the simulations were continued further. In Fig. 1(a)), 
we see only one peak because the state did not decay during our run time. With another 
sequence of random numbers, the meta-stable decays to the ground state. This results in a 
second peak in P(ImL) at ImL ^ as seen in Fig. 1(b). Note that the peak position of 
the meta-stable state is at ImL ~ —0.65 which is same as that in Fig. 1(a). We find that 
P(ReL,ImL) is independent of thermalisation apart from an overall normalisation factor. 
Random processes such as thermalisation can not give rise to such peaks. 

Ideally the system should be allowed to explore all possible states in a single run. If there 
are meta-stable states, the distribution P(ReL,ImL) will have three distinct peaks. The 
ratio of the peak heights of the meta-stable state to the ground state will be given approxi- 
mately by the exponential of the difference in free energies of these two states. Thereby, this 
ratio will be vanishingly small for large volumes. This translates to the vanishing probability 
of observing these meta-stable states starting from any arbitrary initial configuration. These 
states can only be observed by starting from suitable initial configurations. We believe such 
a method of study has no bearing on the properties of the meta-stable states. This is evident 
as we observe almost identical distributions P(ReL, ImL) for different choices of initial con- 
figurations and random number sequences as they all thermalise to a particular meta-stable 
state. 

We have already mention that the physics we study here is different from the meta- 
stability observed near the critical point of a first order transition. For a first order transition, 
meta-stable (stable) state becomes stable (meta-stable) as the system passes through the 
critical point. The meta-stable states are observed close to the transition point. On the other 
hand, the Z(3) meta-stable states exist away from the phase transition point as we discuss 
below. 
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B. Studies with a given volume 



In our simulations for N s = 32, N T = 4, with f3 values up to j3 m ~ 5.80, the system 
thermalises to the ground state irrespective of the initial configurations (listed in the previous 
section) (Fig. 2(a)). The system thermalises to the meta-stable states only for f3 > f3 m 
(Fig. 2(b)). One can argue that the meta-stable states are not observed below f3 m because 
the thermal flucutations are large compared to the barrier height between these states and 
the stable state. In effect this situation can be described by an effective potential without 
any local minima. We mention here that the chiral condensate takes slightly higher value in 
the meta-stable states than that of the ground state. 
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FIG. 2. ReL vs. number of trajectories for three different randomly chosen seeds at (a) /3 = 5.75, 
and (b) (3 = 5.82 for 2-flavor case. 

Higher /3 corresponds to higher temperature and larger quark mass [16, 18]. In the pure 
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FIG. 4. Imaginary vs. real part part of Polyakov loop for two different values of /3. 



577(3) gauge theory, the barrier height between the Z{3) states increases with temperature. 
This should also be true in full QCD, since gluons dominate at higher temperatures. So it 
is expected that the barrier height between the meta-stable state and the absolute ground 
state increases with f3. As a result, the system in meta-stable state should spend on an 
average longer "Monte-Carlo time" in the meta-stable state. Considering different initial 
configurations and random number sequences (for a particular f3 and m g ), we calculate the 
average Monte-Carlo time the system spends in the meta-stable states. We find that this 
time monotonically increases with /3. For example, in Fig. 3 we have shown the case for 
9 = — 27r/3, and f3 = 5.80 and 5.85. For f3 = 5.85, the system stays longer in the meta-stable 
state than f3 = 5.80. For (3 = 5.90, the meta-stable state did not decay during the simulation 
we considered. In Figs. 4, we show ReL vs. ImL for different and f3 values. We see that 
for (3 = 5.90 (Fig. 4(b)), there are no fluctuations between different minima. With various 
different trial configurations and different seeds, we found only two meta-stable states. No 
other (meta)stable state was observed in our simulations, such as 9 = n [7]. 



C. Volume and N T studies 



We extend our studies to different spatial voulmes, V's, and temporal extent, iV T 's, in or- 
der to investigate whether the finite volume effects play any significant role in the observation 
of the meta-stable states. We conside V = 16 3 , 24 3 , 32 3 and N T = 4, 6. For a particular choice 
of (V, N T) rtf and m g ), we always find a /3 m above which there are meta-stable states. For a 
fixed N T . the value of f3 m remains almost the same as volume changed. The average value of 
\L\ in the meta-stable (stable) state remains the same as we increase V from 16 3 — 32 3 as seen 
in Figs. 5. For the different volumes we considered, the Polyakov loop susceptibility does not 
have any dependence on the size of the system. This is expected since at /3 = /3 m , the system 
is away from the transition region. The corresponding correlation length is expected to be 
very small compared to the system size. So we hardly see any finite size effects. In Fig. 5, 
it can be seen that the average value of \L\ and it's fluctuations decrease when A^ r changes 
from 4 to 6. This behavior is observed in both the meta-stable and the ground states. In our 
estimate for /3 m , we find the value increases with iV r . For example, for V = 32 3 we find that 
f3 m ~ 5.80(5.89) for iV r = 4(6). The j3— function [18] results show that a fixed temperature 
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FIG. 5. ReL vs. number of trajectories for V = 16 3 , 32 3 at (a) (3 = 5.85, N r = 4 and (b) /3 = 5.92, 
N T = 6. 



scale corresponds to a higher f3 for larger N T in the /3 range of our study. Interestingly, we 
find that the two /3 m 's obtained for A^ r = 4, 6 correspond to the temperature scales close 
to each other. However, a quantitative comparison between different N T results will require 
higher statistics than the present simulations. 

In all the volumes we study, we find a similar pattern in the behavior of the Polyakov 
loop in the meta-stable and stable states. As expected, f3 m increases with N T . But we would 
like to mention that our study of thermodynamic and continuum limit is far from complete. 
With the increase in the size of the system, the meta-stable states should take longer average 
Monte-Carlo time to decay to the stable state. In our simulations with a small number of 
seeds, we do see that the meta-stable states take longer average MC time to decay with 
the increase in volume. However, a quantitative estimate of the volume dependence of the 
average decay time will need a large number of simulations with different seeds. We mention 
here that same /3 corresponds to lower physical temperature for higher N T . This can lead to 
smaller decay time for the meta-stable states for higher N T (see Fig. 5). 



D. rif and m q dependence 



We have repeated our simulations for 3— flavors of degenerate quarks. In this case, the 
meta-stable states appear at higher values of f3 m compared to that of 2— flavor. Apart from 
higher /3 m , all other findings are similar to those of 2-flavor case. So, we have not shown any 
results for 3— flavor here. We have also considered higher values of (m/T) in our simulations. 
For infinitely heavy quark masses, there is no explicit symmetry breaking of Z(3). So, we 
expect that higher the quark mass smaller is the explicit symmetry breaking. For a larger 
mass am U: d = 0.1, the Monte Carlo histories for /3 = 5.80 has been shown in Fig. 6 together 
with am u ^ = 0.01 for the same /3 and seed. We see that the meta-stable state with am u ^ = 
0.1 decays later than that of am u ^ = 0.01. This suggests that f3 m is smaller for higher (m/T) 
and approaches the pure gauge transition point for m q — >■ oc. 

(m/T) is fixed in our simulations. But for realistic simulations, one should fix the quark 
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FIG. 6. ReL vs. no. trajectories at f3 = 5.80 with masses at am u ^ — 0.01 and am u ^ — 0.1. 

masses. This amounts to considering a smaller (m/T) at higher f3. Smaller (m/T) will 
decrease the barrier height between the meta-stable and ground state. On the other hand, 
the pure gauge effects will make the barrier height grow at higher /3. We expect that the pure 
gauge effects will take over the explicit symmetry breaking term for /3 > /3 m . The value of 
fim we get, should be the lower bound for all the quark masses for which am u ^(f3 m ) < 0.01. 

The quark masses we have considered are larger compared to that of the physical u, d 
quarks. So, if we ignore the strange quark, the f3 m for QCD with d quarks will be higher 
than the value we got from our 2-flavor simulations. Inclusion of strange quark will only 
increase the /3 m . From this, we can conclude that /3 m for realistic case of light u,d quarks, 
and a heavy strange quark will be larger than the value for our 2— flavor calculations. Note 
that f3 m is larger than the /3 C for the Quark-Hadron transition. As a result the meta-stable 
states do not appear near the Quark-Hadron transition but at larger temperatures. 



E. Estimate of temperature scale T m for the meta-stable states 

Now we make an estimate of the temperature T m corresponding to f3 m and am q . This 
requires the /?— function which relates these bare parameters to the lattice cutoff, a. The 
f3— function for Nf = 2 which is appropriate for the action we have used has been studied 
in detail in [18] in the range 5.2 < f3 < 5.6. However, the values of /3 m which is of interest 
to us lie outside this range. But we naively used the same /3— function as in [18] to make an 
estimate of T m . This gives us T m = T(/3 m = 5.80, am q = 0.01) ~ 750 MeV. This temperature 
should be compared with the phase transition temperature which is about 200 MeV. For a 
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better estimate, we should have lattice results for the /3(a) close to ft m . However, we believe 
that this will not drastically reduce T m . As we have argued before, ft m will be larger for 
realistic quark masses with Nf = 2 + 1, 3. So, it is likely that the actual value of T m will be 
larger than our estimate. 

Previous studies have suggested that the explicit symmetry breaking due to quarks is 
small [14]. But as we can see here that the explicit symmetry breaking is strong enough 
to make the Z{3) states with non-zero phase angle unstable below T m . The value of T m > 
750 MeV suggests that this states will not be excited in RHIC. There is a possibility, however, 
that they will be observed at LHC. During thermalization a local region or the whole of the 
fireball may be trapped in one of the meta-stable states. The regions in meta-stable phase 
can then decay through bubble nucleation or via spinodal decomposition depending on 
dynamical evolution of the system [2, 3]. It would be interesting to explore the possible 
consequence of the decay of the meta-stable phases. For temperature above T m the system 
is dominated by gluons and Z{3) is effectively restored. This allows for various types of 
domain walls, static or otherwise [19]. These solutions become unstable below T m . 



V. CONCLUSION 

We have done full QCD simulations with Nf = 2, 3— flavor of degenerate quarks to study 
meta-stable states in the neighborhood of T c . It has been found that the ft value above which 
the meta-stable states appear is close to critical value of ft for the pure gauge confinement- 
deconfinement transition. We estimate the temperature scale to be T m > 750 MeV (for 
Nf = 2) above which meta-stable states can appear. So, we expect that these meta-stable 
states may be observed at LHC. There is a non-trivial change in the shape of the effective 
potential for the Polyakov loop as the system cools through T m . Below T m , there are no 
local minima in the effective potential. So, the meta-stable states, Z 3 domain walls become 
unstable. This may have important consequences for systems such as the fireball in heavy-ion 
collisions or the early Universe. 
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